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ABSTRACT 


Dye  diffusion  experiments  were  made  on  September  11  and  12  1963 

thirty  miles  south  off  Long  Island.  Ten  gallons  of  Rhodamine  B  dye  were 
pumped  at  three  deaths  (0,  3m  and  7m)  in  the  water  each  day.  On  the  first 
day,  the  dye  patch  showed  a  curved  pattern  and  dye  was  diffused  only  above 
10  meters,  perhaps  due  to  rather  calm  weather  conditions.  On  the  second 
day,  the  dye  patch  was  elongated  but  did  not  show  a  remarkable  curvature; 
also,  dye  penetrated  into  deeper  than  70  meters,  owing  to  strong  mixing 
caused  by  winds  of  15  to  20  knots. 

In  Appendix  I,  a  mathematical  model  of  wind-driven  currents  in  the 
ocean  with  a  finite  depth  for  eddy  viscosity  decreasing  exponentially  with 
depth  was  treated.  It  is  shown  that  the  hodograph  of  the  currents  has  a 
stronger  curvature  near  the  surface  for  eddy  viscosity  decreasing  down¬ 
wards  than  the  ordinary  Ekman  spiral  for  a  constant  eddy  viscosity.  In 
Appendix  II,  the  study  of  three-dimensional  structure  of  wind  waves  was 
reviewed.  A  hydrodynamic  model  of  a  stationary  cellular  convection 
superposed  on  a  horizontal  current  was  discussed. 
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I.  INTRODUCTION 


A  series  of  experiments  on  diffusion  of  dye  patches  were  made  by 
a  group  of  oceanographers  of  Lamont  Geological  Observatory  in  various 
parts  of  the  ocean  including  New  York  Bight  and  Bahama  Banks  (Costin 
and  others,  1963),  Among  many  interesting  features  of  dye  patterns 
which  were  observed  in  these  experiments,  curvature  and  striation  are 
the  most  notable.  In  order  to  explain  the  two  features,  a  more  or  less 
consistent  theory  is  yet  to  be  found*  although  the  curved  pattern  is 
attributed  to  the  Ekman  type  spiral  (Gerard  and  Katz,  1963)  and  the 
striation  is  assumed  to  be  caused  either  by  shearing  instability  or  by 
Langmuir  streaks  (Ichiye  and  others,  1964), 

The  experiment  described  here  was  planned  to  get  more  information 
on  mechanism  of  curvature  and  striation  of  dye  patches  by  choosing  an 
experimental  area  off  Long  Island,  as  indicated  in  Fig.  1,  which  is 
relatively  deep  (about  55m),  flat  and  sufficiently  remote  from  the  coast, 
unlike  the  areas  where  the  previous  experiments  were  made.  The 
experimental  procedure  included  releasing  ten  gallons  of  Rhodamine  B 
dye  solution,  tracking  the  dye  patch  with  a  shipborne  Turner  Model  III 
flurometer  and  taking  pictures  from  an  airplane  (Costin  and  others,  1963), 
The  dye  solution  was  pumped  into  the  water  through  three  holes  opened 
at  the  surface,  3m  and  7m  on  a  polyethylene  tube  which  was  suspended 
downward  from  the  boat.  Owing  to  some  difficulties  of  mooring,  no 
current  meter  was  utilized  to  determine  vertical  structure  of  the  current. 
The  surface  drift  of  the  dye  patches  were  measured  by  a  Decca  navigator 
system.  Vertical  distributions  of  temperature  and  salinity  on  the  spot 
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were  measured  several  times  with  a  portable  salinometer.  Winds 
were  measured  frequently  with  a  hand  anemometer  on  the  upper 
deck  of  the  ship  and  also  smoke  signals  were  used  occasionally  to 
indicate  the  direction  of  wind. 

The  research  vessel  used  in  this  experiment  was  R/V  Kyma  of 
New  York  University,  The  ship  proceeded  to  the  experimental  area 
on  the  10th  of  September,  1963  and  the  dye  dumpings  were  made  on 
the  11th  and  12th*  respectively.  The  plane  used  for  aerial  photography 
was  a  Grumman  ivlallarad  on  the  11th  and  a  U.S.  Navy  C-47  on  the  12th. 

II,  Meteorological  and  Hydrographical  Conditions  during  the  Experiment 

The  weather  during  the  10th  and  11th  of  September,  1963  was  rather 
mild.  Typical  weather  charts  of  the  survey  period  are  shown  in  Fig.  2. 
The  cold  front  passed  the  New  York  area  at  about  0500  of  September  10. 
The  state  of  the  sea  over  the  experimental  area  was  calm  from  noon  to 
about  1600  and  swells  of  1/2  meter  high  and  50  meters  long  came  from 
the  southeast.  At  about  1630  the  southeast  wind  reached  2  Beaufort  scale. 
On  September  11th  the  area  was  covered  by  a  high  pressure  area  and 
generally  the  winds  were  moderate.  The  wind  speed  and  direction 
measured  on  board  the  ship  are  listed  in  Table  1.  The  directions  of 
wind  waves  were  almost  the  same  as  the  wind  and  their  heights  were 
about  0.7  meters  throughout  the  day.  However,  the  swells  coming  from 
the  south  superposed  on  the  wind  waves  and  thus  caused  cross  seas  after 
the  shift  of  wind  direction  from  east  to  south. 

During  the  night  of  September  11th  the  southerly  wind  increased  its 
speed.  During  the  day  of  the  12th  the  wind  reached  25  to  40  knots.  Both 
wind  waves  and  swells  came  from  south  and  their  height  reached  about 


‘ 


:1 


•  ' 


.  . 


i.x, 


* 


v  >  :i 


i 


>  *  •  •  ; 


* 


*■ 


i 


.  * 


r.>  • 


? 


-3- 


four  meters  in  the  afternoon.  These  strong  southerly  winds  were  caused  by 

a  cold  front  approaching  from  the  west. 

Table  1  Wind  Measurements 


Date_ 
Sept.  11 


Sept.  12 


T  im  e 

Wind  Speed  (kt. ) 

Wind  Direction 

0845 

18 

E 

1140 

12 

ESE 

1325 

12 

E 

1415 

11 

ESE 

1445 

14 

SE 

1645 

12 

S 

1750 

10 

s 

1855 

12 

SSE 

0620 

1 6  — - 1 9 

SSE 

0945 

18  v>21 

SSE 

1040 

20  -  2  3 

SSE 

11  30 

22  .-v  25 

S 

1155 

23  r-j  2 6 

5 

1415 

35  v,40 

S 

Temperature  and  salinity  profiles  obtained  with  a  portable  salinometer 
are  plotted  in  Fig.  3.  It  is  noted  that  the  thermocline  existed  below  25 
meters  on  both  days  but  temperature  jump  at  the  thermocline  diminished 
conspicuously  on  the  12th  owing  to  warming  of  deep  water  due  to  the  wind 
stirring.  Salinity  distributions  indicate  that  there  was  the  halocline  at  the 
same  depth  as  the  thermocline  on  the  11th  but  it  disappeared  on  the  12th. 
III.  Patterns  of  Dye  Patches 

Aerial  photographs  of  dye  patches  were  taken  at  several  minutes' 
interval  with  a  polaroid  or  a  Graflex  camera  using  Kodak  Wratten  23A 
filters.  Most  pictures  contain  not  only  a  dye  patch  but  also  the  ship,  so 
that  the  dimension  of  dye  patch  may  be  estimated.  Also,  the  heading  of 
the  ship  was  recorded  almost  every  minute  and  the  orientation  of  the 
patch  can  be  determined  with  reasonable  accuracy.  Several  pictures 
contain  the  smoke  signal  released  aboard  the  ship  and  the  relationships 
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between  wind  and  dye  patch  orientation  were  determined.  Sketches 
of  the  dye  patches  based  on  aerial  photographs  are  shown  in  Fig.  4. 

On  the  11th,  dye  was  dumped  at  1000  E.S.T.  This  patch  is  called  the 
first  patch  hereafter.  Striations  became  visible  about  twenty  minutes 
after  the  dumping.  The  wind  direction  was  east  cr  east-southeast 
during  the  time  of  aerial  reconaissance .  The  striations  were  almost 
parallel  to  the  wind  direction,  along  which  the  furrows  of  short-crested 
waves  were  aligned.  At  1040  the  elongation  of  the  patch  became  obvious 
and  the  striation  became  more  regular.  At  1050  the  patch  showed  a 
clockwise  curvature  from  head  to  tail.  (At. 'the  head,  concentration  is 
high,  while  at  the  tail,  the  concentration  is  low).  After  11:30  the 
patterns  of  the  patch  became  almost  stable  although  the  size  increased 
as  the  time  proceeded.  The  curvature  became  striking.  The  striations, 
were  in  the  east-west  direction,  almost  parallel  to  the  wind  and  had 
rather  regular  spacings  of  50  to  100  meters.  Until  about  1200,  the 
head  was  broad  and  the  tail  was  narrow  but  after  then,  the  tail  became 
broader  than  head  and  the  size  of  the  patch  increased  rapidly. 

On  the  12th  the  dye  was  dumped  at  0940  E.S.T.  This  patch  is 
called  the  second  patch.  Striations  began  to  form  almost  ten  minutes 
after  the  dumping,  earlier  than  on  the  previous  day.  At  1015  the 
striations  became  very  conspicuous  and  their  orientation  was  north- 
south,  with  spacings  of  20  to  50  meters.  This  orientation  was  also 
almost  parallel  to  the  prevailing  wind.  At  the  same  time,  the  whole 
dye  patch  began  to  elongate  in  the  direction  slightly  oblique  to  the  wind. 
This  elongation  became  more  and  more  conspicuous  and  its  direction 
was  deviated  to  the  right  of  the  wind,  with  angles  between  20  and  30 
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degrees.  At  about  1050  larger  striations  with  spacings  of  100  meters 

<■  i  '  •'  t  [  , 

or  more  began  to  develop.  After  1130  the  pattern  of  the  patch  became 

,/  ! 

rather  stable  except  some  development  of  larger  striations.  The  size 
also  gradually  increased  but  the  increase  was  not  so  conspicuous  as  on 
the  previous  day. 

IV  Vertical  profiles  of  concentration  and  drift  of  the  patches 

The  profiles  of  concentration  of  dye  measured  with  the  fluorometer 
are  plotted  in  Fig,  5.  Since  no  buoy  system  was  available,  such  measure¬ 
ments  were  made  only  at  head,  middle  and  tail  of  the  dye  patches.  This 
figure  clearly  indicates  that  the  dye  on  the  11th  spread  only  above  10  m 
but  on  the  12th  23  m  or  deeper.  Also,  in  the  first  patch  the  dye  reached 
almost  the  same  depth  from  head  to  tail  (axis  determined  from  the  first 
moment  of  the  concentration  about  the  depth  is  3,0m,  3.  lm  and  3*7m  at 
the  head,  the  middle  part  and  the  tail  respectively),  but  in  the  second 
patch  the  head  seemed  to  reach  deeper  than  the  tail.  Location  of  dye 
patches  determined  with  a  Decca  navigation  system  is  plotted  in  Fig.  1. 

On  the  second  day,  the  drift  of  the  dye  patch  was  to  the  southeast  be¬ 
tween  0900  and  1300  opposite  the  direction  of  the  wind.  This  movement 
seems  to  be  due  to  the  tidal  current.  The  tides  *t  Shinnecock  Inlet 
(40°  50‘  N,  72<*28l,  south  side  of  Long  Island)  on  September  12  is  as 
follows:  HW,  0137;  LW,  0800;  HW,  1427;  LW,  2053.  The  southward 
flow  seems  to  correspond  to  the  flood  tide.  There  is  no  data  of  tidal 
currents  in  this  area  and  it  is  hard  to  confirm  the  relationships  between 
the  tide  and  tidal  currents.  However,  since  the  tidal  currents  are 
almost  uniform  vertically  except  close  to  the  bottom,  their  effe  ct 
is  to  transport  the  entire  part  of  the  dye  patch. 
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Therefjre,  in  the  present  experiment  dealing  with  a  moving  source,  the 
pattern  of  the  dye  patch  might  not  be  disturbed  by  presence  of  tidal 
currents. 

In  Fig.  4  the  directions  of  smoke  and  wind  direction  and  speed 
measured  with  a  hand  anemometer  were  inserted  whenever  such  data 
were  available.  Also,  the  time  indicated  with  parenthesis  means  that 
a  photograph  at  that  time  did  not  show  the  ship  and  thus  the  size  of  the 
patch  was  estimated  by  interpolation. 

Photograph  1  shows  the  first  patch  at  1255  on  September  11th  and 
Photograph  2  shows  the  second  patch  at  1240  on  the  12th.  The  altitude 
at  which  the  pictures  were  taken  was  1300  m  for  the  first  patch  and 
650  m  for  the  second  patch, 

V  Theoretical  Discussions 

One  of  the  striking  features  cf  the  dye  patches  at  this  experiment 
is  that  the  first  patch  showed  a  conspicuous  curved  tail,  while  the  second 
patch  elongated  but  was  rather  straight.  The  vertical  distributions  of  dye 
concentration  in  Fig.  5  indicate  that  dye  was  found  only  above  10  meters 
in  the  first  patch  but  it  penetrated  deeper  than  20  m  in  the  second  patch 
suggesting  that  the  vertical  mixing  was  more  intense  on  the  second  day 
owing  to  the  stirring  effect  of  winds.  The  warming  of  water  below  the 
thermocline  on  the  second  day  (Fig.  3)  indicates  that  the  effect  cf  stir¬ 
ring  reached  below  the  thermocline. 

The  lack  of  curvature  in  the  second  patch  may  be  due  to  the  increase 
of  the  vertical  eddy  viscosity  by  the  effect  of  wind  stirring.  The  depth 
of  frictional  resistance  D  for  a  constant  eddy  viscosity  equals  to  [2  ^  ) 

where  /  •  is  the  eddy  viscosity  and  X-  is  the  Coriolis'  parameter. 
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(Sverdrup  et  al,  The  Oceans,  pp,  492-493).  When  the  eddy  viscosity 
increases,  D  also  increases.  If  D  becomes  larger  than  the  depth  of 
the  sea,  the  Ekman  spiral  does  not  develop  and  instead  the  velocity 
vector  tends  to  have  almost  the  same  direction  throughout  the  depths, 
although  the  directions  are  reversd  near  the  bottom  (Defant,  1961, 
p.  402).  This  prediction  based  on  the  classical  Ekman  theory  of  the 
wind  current  in  the  sea  of  a  finite  depth  can  be  applied  to  the  second 
patch.  It  was  elongated  to  the  direction  deviated  clockwise  by  about 
20  to  30  degrees  from  the  wind  and  there  was  no  curvature. 

The  curvature  of  the  first  dye  patch,  which  was  clearly  recognized 
almost  one  hour  after  the  dumping,  seemed  to  be  due  to  the  Ekman 
spiral,  when  the  eddy  viscosity  and  thus  the  friction  depth  D  are 
sufficiently  small,  so  that  the  effect  of  the  bottom  did  not  influence 
the  structure  of  the  wind  current.  In  fact,  the  wind  of  8  to  12  knots 
was  almost  constantly  blowing  from  east  southeast  or  east  since  the 
afternoon  of  the  10th  of  September.  Owing  to  the  moderate  wind,  the 
eddy  viscosity  was  not  so  large  as  in  the  second  patch  and  thus  the 
Ekman  type  wind  currents  were  well  established  during  the  morning 
of  the  11th.  This  is  the  reason  why  the  curvature  of  the  dye  patch  so 
quickly  developed. 

When  the  eddy  viscosity  decreased  rapidly  with  depth,  the 
curvature  of  the  hodograph  of  the  Ekman  spiral  becomes  larger  than 
for  the  uniform  eddy  viscosity,  particularly  near  the  surface  (Fjeldstat, 
1929).  The  eddy  diffusivity  decreasing  with  depth  also  makes  vertical 
diffusion  of  dye  concentrate  near  the  surface  as  manifested  by  the 
first  patch.  The  analytical  discussion  on  the  drift  current  due  to  the 
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eddy  viscosity  variable  with  depth  is  presented  in  Appendix  1. 

The  patterns  of  dye  patches  described  in  Fig.  4  indicate  conspicuous 
striations  which  were  almost  parallel  to  the  wind  direction  on  the  second 
day.  On  the  first  day,  the  striations  were  parallel  to  the  wind  which  was 
easterly  at  first,  but  after  the  wind  shifted  to  the  south,  reorientation  of  tine  striatiom 
did  not  immediately  follow. 

The  orientation  of  striations  suggest  that  it  might  be  caused  by 
surface  streaks  which  were  explained  as  the  effect  of  surface  film  by 
Welander  (1962).  However,  the  wide  spacings  (a  few  hundred  meters) 
and  slow  response  of  alignment  to  wind  directions  are  quite  contrary  to 
the  conditions  envisaged  by  7/elander. 

The  process  of  development  of  striations  will  be  described  in  the 
following,  according  to  visual  observations  on  board  the  ship  and 
inspection  of  aerial  photographs.  On  September  11th,  wind  waves  came 
from  east  in  the  morning  due  to  the  easterly  wind  as  described  in 
Section  2.  Aerial  photographs  at  1022,  1040  and  1127  indicate  that  the 
furrows  or  streaks  with  a  width  of  about  a  ship's  length  (22  meters) 
were  superposed  on  this  wave  pattern,  ran  from  east  to  west  and  were 
orthogonal  to  the  crest  lines  of  waves  which  have  crest  lengths  of  three 
to  four  ships'  length.  No  definite  meteorological  disturbance  which 
might  generate  such  swells  is  seen  in  weather  charts  of  Fig.  2,  Also, 
these  furrows  or  streaks  appeared  to  be  almost  stationary.  Close 
inspection  of  the  photographs  mentioned  above  shows  that  the  furrows 
seem  to  correspond  to  the  dark  bands  between  the  striations.  This 
suggests  that  the  furrows  were  associated  with  the  divergence  or  con¬ 
vergence  of  the  surface  current. 
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On  the  12th,  striation  began  to  develop  almost  20  minutes  after  the 

i  1  *  ’  '  *  ' 

dumping.  Although  the  furrows  parallel  to  the  wind  were  recognized, 

jj  i  1  , 

they  were  not  so  clear  as  on  the  previous  day.  The  spacings  of  striations 

i  ,  |  )  i 

were  irregular  and  smaller  than  a  ship’s  length  at  first.  The  striations 

l 

were  developed  leeward  in  the  subsurface  part  of  the  dye  patch  as  seen 
from  the  direction  of  smoke  contrary  to  the  case  of  the  previous  day. 
Fingers  of  the  dye  patch  seemed  to  be  elongated  in  the  direction  of 
currents  produced  by  breaking  of  wave  crests.  After  40  to  50  minutes 
of  the  dumping,  the  fingers  developed  windward  and  spacings  reached 
several  ships’  length.  In  this  stage,  the  direction  of  elongation  of  the 
whole  patch  was  deviated  to  the  right  of  wind  direction  as  stated  in  Section 
3,  but  the  striations  were  oriented  almost  parallel  to  the  wind. 

It  is  quite  certain  that  the  three  dimensional  structure  of  waves  or 
short-crested  waves  are  important  to  the  generation  of  striations.  As 
Phillips  (1957)  and  Lighthill  (1962)  discussed,  turbulence  of  wind  may  be 
an  important  factor  to  wind  waves  in  the  developing  stage.  As  explained 
in  Appendix  II,  shorter  waves  including  capillary  wave  regime  propagate 
statistically  in  two  main  directions  with  angles  more  than  60°  to  the  mean 
wind  direction.  The  interference  among  these  wave  trains  seems  to 
produce  alternatively  rough  surfaces  along  the  crests  whose  width  is 
more  than  several  wave  lengths  and  smooth  surfaces  along  the  furrows 
which  are  aligned  with  the  mean  wind  directions.  The  differential 
roughness  along  the  (short)  crests  and  along  the  furrows  is  twofold. 

The  large  roughness  along  the  crests  causes  increase  of  the  shear  in 
the  surface  layer  of  the  water  and  thus  faster  surface  currents  than 
along  the  furrows,  occasionally  reinforced  by  breakings  of  tips  of 
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waves.  This  produces  fingers  of  the  dye  patch.  Also,  the  reduced 

/ 

I  ;  I  , 

shear  along  the  furrows  due  to  the  smooth  surface  causes  the  diver - 

i 

gence  or  convergence  of  the  surface  current  there  (Appendix  II), 

;  I  ,  •  •  i  i  = 

Along  the  divergence  line,  the  water  with  less  concentration  of  dye 
is  brought  from  below  to  the  surface,  thus  forming  the  dark  bands 
between  the  striations. 
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Appendix  I  -  W  ind-driven  Current  for  an  Eddy  Viscosity  Variable  with  Depth 

Effects  of  surface  waves  on  generation  of  turbulence  in  the  water 
were  discussed  semi-empirically  by  Dobronklonsky  (1947)  and 
Shebalin  ( 1957)Ulchiye  (1962).  The  heuristic  result  of  these  studies  is 
an  expression  of  the  eddy  viscosity  due  to  waves.  When  the  depth  is 
finite,  the  eddy  viscosity  is  given  by 


1  *  ii  iv!  ('  ^  ^ ) _ _ 

iv  K-  ii  u.$h2^fc(hi  -z) 


(i) 


where  kQ  is  Karman's  constant  (=0.4),  p  the  density  of  the  water,  /:- 
the  wave  height,  T  the  period,  k  is  the  wave  number,  z  is  the  depth 
below  mean  sea  level  and  H  is  the  depth  of  the  water.  When  the  water 
is  deep  compared  with  the  wave  length,  this  equation  becomes 


These  two  relations  indicate  that  in  general,  the  eddy  viscosity 
increases  with  increasing  height  and  frequency  and  that  the  ratio  of  its 
decrement  with  depth  becomes  larger  as  increasing  wave  number. 

In  order  to  examine  dependence  of  the  Ekman  wind  current  on  differ¬ 
ences  of  the  eddy  viscosity,  the  simple  case  of  steady  state  and  pure  drift 
current  is  treated.  Also,  for  mathematical  simplicity,  the  eddy  viscosity 
of  the  type  of  equation  (2)  is  assumed,  though  the  depth  is  taken  finite. 

The  equations  of  motion  are  given  by 


(3) 


where  W  =  u  +t  u  and  v  being  x-  and  y-  components  of  current, 

respectively.  The  boundary  conditions  at  the  surface  (z  =  0)  and  at  the 
bottom  (z  =  h)  are  expressed  by 


(4) 
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W  =  (5) 

in  which  the  wind  stress  is  assumed  to  have  y-  component  only. 

When  the  variable  (z)  is  changed  into  x,  using  the  relationship 

>1-  (6) 


equation  (3)  becomes 
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>c  dL  ( - x  7-  C{  !/1/ 
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d  X 


,  6 <  ~  ~ 
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iv 


(7) 


where  •/ -X  _  f  ,n  yA\"/.  The  boundary  condition  (4)  is  rewritten  as 

r  -  ;  i'/b-A 


d  w 

c  I  X 


<•  ") 
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-  / 


(8) 


where  T  =  /X  n  \~/.  The  solution  of  equation  (6)  can  be  expressed  by 

-  /  r  j 

cylinder  functions  like 

J-  '  A  .  ' f  >  0  .*L\  ,  n  i  i1^  / 

•  t  AM,  ( A *  -OtBH  i  r  -)  ■;  (9) 
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j  i  ^ 


where  p  =  (  )and  H  (  and  rT  are  Hankel’s  function  of  the 

first  and  the  second  class,  respectively.  Integration  constants  A  and  B 
can  be  determined  from  the  boundary  conditions  at  the  surface  (x  =  1)  and 
at  the  bottom  (x  =  a  =  exp  The  conditions  (8)  and  (5)  respectively 

lead  to  the  following  relationships: 
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(ii) 


Substitution  of  the  values  of  A  and  B  determined  from  (10)  and  (11)  and 
(9)  yields 
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where 


n 


-» 


» 


V 


(13a.b) 


When  ~i  j'j  is  much  less  than  unity,  the  eddy  viscosity  is  almost 
uniform  with  depth  and  the  solution  for  wind-driven  currents  is  al¬ 
most  the  same  as  the  original  one  by  Ekman.  Therefore,  the  inter¬ 
esting  case  is  that  the  eddy  viscosity  decreases  rapidly  with  depth 
or  the  parameter  | )  is  much  larger  than  unity.  When  the  depth 
is  taken  as  50  m,  equation  (2)  yields  =  25  (cm^/sec)  _  _  jq 

for  waves  with  wave  length  13.  6  m,  wave  height  0.  5  m  and  period  3 
sec.  With  these  values  of  Yj  .  and  ^  we  have 


(14) 


Therefore,  m.  can  be  expressed  by  an  expansion  of  Hankel  function  for 
a  small  value  of  argument  while  m  is  given  by  its  asymptotic  expansion. 


When  only  the  first  order  term  is  retained,  we  have 

"n)j  -  <  (it)  hfs  ,  (j  =  H ) 


(15) 
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where  r  9  p  ^  ~~  h,  and  *s  Euler's  constant  (  =  1.781)  (Jahnke  and 
Ende,  1945),  Equations  (15a)  and  (15b)  indicate  that  j  t  (  j  j  y^  ..  j  . 
The  approximate  expression  for  W  in  the  depth  range  satisfying  o(g<l 


is  given  by 
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It  is  noted  that  the  expression  (17)  does  not  depend  on  z  explicitely. 
The  value  of  cor  responding  to  p  =  0.04  given  by  (14)  equals  0.21 
or  about  12°.  Therefore,  the  surface  drift  deviates  from  the  wind 
to  the  right  by  about  78°,  in  contrast  to  the  classical  Ekman  drift. 

When  the  depth  satisfies  the  relation  ,  the  approxi¬ 

mate  expression  for  W  becomes 


N 


-  r  A 
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where  r  =  £  p  ^  "E  *  .  The  distribution  of  current  vectors  with 

depth  for  the  same  numerical  values  of  //  •  i  and  j:,  as  leading 
to  (14)  are  plotted  in  Fig.  6.  It  is  remarkable  that  the  decrement  of 
velocity  with  depth  is  quite  rapid  and  that  the  curvature  of  the  resultant 
spiral  is  quite  conspicuous  near  the  surface.  For  comparison,  the 
vertical  structure  of  the  wind-driven  current  for  the  constant  eddy 
viscosity  is  shown  in  the  curves  taken  from  Defant  (1961)  in  Fig.  5. 

In  these  curves,  D  =  <7"  j  t^ie  frictional  depth  and  thus  D 

equaling  to  4h,  2h  or  0.8h  yields  the  eddy  viscosity  of  2  x  10^,  5  x  10^ 
or  80  (cm^/sec),  respectively,  when  h  =  50  m.  All  these  curves 
indicate  that  the  curvatures  of  the  spiral  arenot  so  strong  near  the 
surface  as  in  the  previous  example. 

Kozlov  (1963)  treated  the  similar  problem  with  eddy  viscosity 
varying  with  depth  according  to  a  power  law  like  ( !  i  -  z )  1 

His  result  (Fig.  1  in  his  paper)  also  indicates  that  the  curvature  of 
the  spiral  becomes  more  conspicuous  as  n  increases,  although  his 
curves  show  larger  curvature  towards  the  bottom. 
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Appendix  II  ^Mechanism  of  Striations  of  Dye  Patches 

At  the  present  time  there  are  two  kinds  of  theory  explaining  generation 
of  waves  by  wind.  One  is  a  model  considered  by  Phillips  (1957)  who  dis¬ 
cussed  the  random  pressure  fluctuations  associated  with  turbulent  wind 
as  a  main  cause  of  wind  waves.  His  theory  indicated  that  waves  with 
length  /\  would  be  resonantly  generated  if  an  angle  (-)  between  their 
direction  of  propagation  and  that  of  convectional  pressure  pattern  having  a  speed  u 
satisfied 


where  C  is  the  phase  velocity  (  =  ju  )• 

Another  theory  developed  by  /files  (1957)  explained  the  growth  of 
waves  as  an  instability  in  the  shearing  flow  of  air -water  system.  This 
model  is  based  on  the  purely  two-dimensional  structure  of  the  mean  wind 
which  varies  only  in  the  vertical  direction.  The  turbulence  of  the  wind 
plays  no  essential  part  in  this  model.  Later  Miles  (I960)  combined  these 
two  models.  He  showed  that  in  the  initial  stage  of  wave  generation  the 
random  pressure  excited  the  interface  oscillations  which  are  then  augmented 
by  the  instability  of  shear  flow.  Longuet-Higgins  (1963)  showed  that  Miles' 
model  is  more  realistic,  since  a  turbulent  component  in  the  mean  square 
air  pressure  was  at  least  two  orders  of  magnitude  smaller  than  postulated 
in  Phillips'  theory.  From  measurements  of  vertical  displacement  and 
slopes  of  sea  level  with  a  floating  buoy,  he  determined  the  directional 
spectrum  /_  /o 

defined  so  that  j-  f  '  (j'\  (jC(j '  *s  the  contribution  to  the  mean  square 
value  of  surface  displacements  due  to  wave  elements  which  have  ranges  of 
frequency  (  '' f  '/■  t <//.'  )  and  direction  (  (j  .  cj)  i  <!  (j>  )  . 
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of  the  surface  waves.  This  directional  spectrum  is 
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He  then  determined  an  approximate  r.m.s.  angle  'C; 
direction  of  wave  propagation  and  that  of  the  wind  for  different 
values  of  6 "  .  The  values  of  ,7/  obtained  from  a  measured 

spectrum  is  in  good  agreement  with  files'  formula  COS  ^(c  / Jj. . _  j » 

where  c  is  the  wave  velocity  (  =  J  ']  'A /2  IT  )  an<*  ~0~'U*  *s 

(7 

wind  velocity  at  a  height  of  0.2  wave  lengths.  Thus,  for  the  waves 


with  periods  shorter  than  3.3  seconds,  the  values  of  are  within 
a  range  of  50°  to  70°,  slightly  increasing  with  decreasing  periods, 
while  for  the  longer  waves,  the  values  of  'd  decreases  rapidly 


with  increasing  periods  or  wave  lengths.  Therefore,  the  shorter 
waves  which  contribute  to  hydrodynamic  roughness  of  the  sea  sur¬ 
face  have  "two-peaked"  directional  spectrum.  This  may  suggest 
that  the  two-dimensional  structure  of  sea  surface,  such  as  alternative 
smooth  and  rough  portions,  is  attributed  to  the  spectral  character¬ 
istics  of  the  wind  waves. 

A  simple  mathematical  model  will  be  developed  for  cell-like 
motion  induced  by  wind  blowing  over  alternative  smooth  and  rough  sea 
surface.  Since  the  roughness  of  the  sea  surface  changes  wind  stresses, 
the  effect  of  variable  roughness  causes  variable  shearing  motion  of  the 
water  near  the  surface.  It  is  assumed  that  the  wind  blows  in  x-direction 


i  which  is  constant  in  the  x-direction.  The 
z-axis  is  taken  positive  downwards,  the  y-axis  is  perpendicular  to  the 


causing  a  shear  flow  w  (  J  y  ) 

\  J  >  "  / 


x-axis  and  the  origin  is  taken  at  the  sea  surface.  Only  the  mean  motion 


is  considered.  Then  equations  of  motion  are 
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where  v  and  w  are  velocity  components  in  the  y-  and  z-  direction 
respectively  and  J/  is  the  kinetic  viscosity.  Since  it  is  assumed  that 
the  motion  is  uniform  in  the  x-direction,  all  quantities  are  independent 
on  x.  The  equation  of  continuity  becomes 
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(4) 


This  equation  yields  a  stream  function  for  v  and  w,  such  as, 
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The  equation  for  vorticity  K  (  =  n  ^  T-  )  obtained  by  cross- 

differentiation  of  (2)  and  (3)  is 
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and 
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All  quantitities  in  equation  (6)  are  dimensionless  and  R  is  Reynolds' 


number  defined  as  (\  j  [J  ,  where  a  and  vQ  are  characteristic 


length  and  velocity,  respectively.  If  we  consider  a  equaling  100  m  as 


the  spacing  distance  of  the  striations  and  assume 


r 


to  have  an 


order  of  1  cm/  sec,  the  order  of  R  =  10  for  the  molecular  viscosity 
equaling  0.  1  cm^/ sec.  Therefore,  equation  (6)  and  corresponding 
equation  (1)  can  be  expanded  in  a  series  of  powers  of  R“*. 


(Since  the  wind-driven  current  f  }  is  turbulent,  it  might  be  con¬ 


tended  that  eddy  viscosity  should  be  used  for  J  -1  instead  of  molecular 
viscosity.  Even  so,  the  eddy  viscosity  has  an  order  of  magnitude  of 
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2  2/  j 

10  cm  /sec  when  Thorade's  formula  (Defant;  1961,  p.  104)  for 
eddy  viscosity  due  to  wind  current  is  applied,  or  of  10  cm  /  sec 
when  Richardson's  4/ 3-power  law  (Defant,  ibid,  p.  390)  is  applied. 
Therefore,  R  has  an  order  of  10  to  10^  and  the  expansion  with  R“* 
seems  to  be  reasonable.) 

The  zeroth  order  term  of  equation  (10)  is 
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(8) 


Corresponding  term  for  equation  (1)  is  obtained  by  substituting 


U 


into  of  (8).  This  indicates  that  the  velocity  u  and  vorticity 

are  expressed  by 


-■  'vjy 


h  "V  ) 


(9) 


where  f  and  are  arbitrary  functions. 

When  a  rectangular  cell  is  considered,  the  boundary  conditions  for  ijy 

I 

are  as  follows: 


'I  N 


-  -  ^  r/-~  *. 


z 


at  y  =  0  and  y  =  1 


(10) 


at  z  =  0  and  Z  -  d  (11) 


where  d  is  a  dimensionless  depth  of  the  lower  boundary.  These  conditions 


can  be  represented  by  a  simple  relation 


at  y  =  0  and  y  =  1 
and 

z  =  0  7  =  d  (12) 


Immediate  result  from  equation  (9)  and  (12)  is  that  at  the  surface  (  z  =  0) 


In  other  words,  the  vertical  shear  of  /  j  at  the  surface  is  proportional  to 
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the  lateral  velocity  at  the  surface.  This  indicates  that  the  vertical 
shear  is  the  strongest  at  the  center  of  a  cell  with  one  gyre  and  it 
vanishes  at  the  boundaries  y  =  0  and  1,  (Fig.  7)  Therefore,  the  slicks 
which  may  reduce  shear  near  the  surface  may  correspond  to  either 
divergence  or  convergence. 

As  a  simple  example,  the  function  JT(  <j,  )  is  assumed  as  a  linear 
function  of  U^,  Then  the  equation  for  x  of  (9)  becomes 


"V — 7  2 


V  ‘  y  -  a  -\y  -  b  <14> 

By  choosing  a  and  b  arbitrarily,  various  types  of  streamlines  can  be 
obtained  from  the  solution  of  (14)  under  the  boundary  condition  (2).  The 


t  s  * 


solution  is  expressed  by 
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In  expression  (15),  the  sine  terms  including  should  be  replaced  with 


terms  such  as  _v  ,  f,  x  -  ^  etc,  if  a  is  negative.  Some  examples  of 


streamlines  determined  from  equation  (15)  are  shown  in  Fig.  7. 
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